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ENDOCARDIAL MAPPING SYSTEM In general the invention provides a system and method 

which permits the location of catheter electrodes to be 

TECHNICAL FIELD visualized in the three-dimensional map. 

Hie invention discloses the apparatus and technique for < The invention may also be used to provide a two- 
forming a three-dimensional electrical map of the interior of 5 dimensional map of electrical potential at or below the 
a heart chamber, and a related technique for forming a myocardial tissue surface. 

two-dimensional subsurface map at a particular location in Additional features of the invention will appear from the 

the endocardial wall. following description in which the illustrative embodiment 

is set forth in detail in conjunction with the accompanying 

BACKGROUND ART 30 drawings. It should be understood that many modifications 

to the invention, and in particular to the preferred embodi- 

It is common to measure the electrical potentials present ment illustrated in these drawings, may be made without 

on the interior surface of the heart as a part of an electro- departing from the scope of the invention, 

physiologic study of a patient's heart. Typically such mea- r,.^ . . . . c 
r J f • }• # j • i e ir FIG. 1 is a schematic view of the system, 

surements are used to form a two-dimensional map of the 15 ' 

electrical activity of the heart muscle. An electrophysiologist FIG 2 15 a view of the catheter assembly placed in an 

will use the map to locate centers of ectopic electrical endocardial cavity. 

activity occurring within the cardiac tissues. One traditional FIG. 3 is a schematic view of the catheter assembly. 

mapping technique involves a sequence of electrical mea- FIG. 4 is a view of the mapping catheter with the 

surements taken from mobile electrodes inserted into the 20 deformable lead body in the collapsed position. 

heart chamber and placed in contact with the surface of the FIG 5 ^ a view of the ma p ping cat heter with the 

heart. An alternative mapping technique takes essentially deformable lead body in the expanded position. 

simultaneous measurements from a floating electrode array rrp ^ . c c ... 

,. , r , . • i . , FIG. 6 is a view of the reference catheter, 
to generate a two-dimensional map of electrical potentials. 

~, ... . r*t_ i ■ i * i * m FIG. 7 is a schematic view representing the display of the 

The two-dimensional maps of the electrical potentials at /D three-dimensional ma 

the endocardial surface generated by these traditional pro- ree 111160510111 ma P- 

cesses suffer many defects. Traditional systems have been FIG 8 155 a side view of an alternate reference catheter, 

limited in resolution by the number of electrodes used. The FIG. 9 is a side view of an alternate reference catheter, 

number of electrodes dictated the number of points for FIG. 10 is a perspective view of an alternate distal tip. 

which the electrical activity of the endocardial surface could 30 FIG. 11 is a schematic view representing the display of the 

be mapped. Therefore, progress in endocardial mapping has subsurface two-dimensional map. 

involved either the introduction of progressively more elec- FIG n & a schematic flow chart of the steps in the 

trodes on the mapping catheter or improved flexibility for method, 
moving a small mapping probe with electrodes from place to 

place on the endocardial surface. Direct contact with elec- 35 DETAILED DISCLOSURE 

trically active tissue is required by most systems in the prior T , 4 , t Ci , , . , f 

. « A .7 ,. / i * In general, the system of the present invention is used for 

art in order to obtain well conditioned electrical signals. An ^ ^ ac(ivi of ^ surface of & 

exception is a non-contact approach with spot electrodes hear( chamber 8Q The ma . catheter assemb , u 

Ttee spot electrodes spatially average the electrical signal a flexib , e , ead b n ^ a deformable 

through their conical v.ew of the blood media This «° dislal lead body 74 . ^ deformab i e distal lead body 74 caD 

approach therefore also produces one signal for each elec- . c j ■ * * ui u 

* rn „ . /. . , u . j j. , be formed into a stable space filling geometric shape after 

trode. The small number of signals from the endocardial • * j *• ■ ♦ .u i_ * % on ^fL- j c i_i j- * i 

„ .„ «. • •»..,•* . .i i .i_ mtroduction mto the heart cavity 80. This deformable distal 

wall will result in the inability to accurately resolve the 1 , . , - A • , , i * j m j c • 

. - 7 _ - t lead body 74 includes an electrode array 19 defining a 

location of ectopic issue masses. In the pnor art, iso- number of electrode sites . Th e mapping calheter assemb i y 

potentials are interpolated and plotted on a rechlinear map „ alsQ indudes a reference e]ectrode ferab] laced on 

which can only crudely represent the unfolded intenor a reference u M fa . , umen 

surface of the heart. Such two-dimensional maps are gen- g2 formed ^ ^ &exMe , ead bod ?2 & d ^ ^ lead 

erated by mterpolation processes which fill in contours , , - A ~- c 4l _ . li i*u j- * i *■ 

, , 3 • . . • p i.j body 74. The reference catheter assembly 16 has a distal tip 

based upon a limited set of measurements. Such interpolated . \ . ~ A . , V . . \ 

j. . , . . n A , - . . -r <?n electrode assembly 24 which may be used to probe the heart 

two-dimensional maps have significant deficiencies. First, if 3U , , J , . . . \ L1 

^ -la. . 1 j • wall. This distal contact electrode assembly 24 provides a 

a localized ectopic iocus is between two electrode views c , . . , f r i u* tt_ u i 

, ,„ r - r , - surface electrical reference for calibration. The physical 

such a map will at best show the ectopic focus overlaying , h of , he reference ca , he(er M (aken ^ tfae iUon of 

both electrodes and all points m between and at worst wdl (he e , ectrode affa 19 hef ide , reference whjch 

not see it at all. Second, the two dimenstonal map, since it fae ^ ^ e „ ^ e refer . 

contains no chamber geometry information, cannot indicate 5:> 3 #u . 1/: , , , .,. r*u i ♦ ^ 

• i u • *u *i_ i • . . , f4l _ . ence catheter 16 also stabilizes the position of the electrode 

precisely where m the three dimensional volume of the heat 1A , . , . U1 

r . . ' , . . . . , . , t , m. . 4 array 19 which is desirable. 

chamber an electrical signal is located. The inability to ' . , . , . , 

accurately characterize the size and location of ectopic tissue ™ e *f str ^ ctural elements provide a mapping catheter 

frustrates the delivery of certain therapies such as"ablation". ass f mbl y which be readily positioned within the heart 

60 and used to acquire highly accurate information concerning 

SUMMARY DISCLOSURE ^ e e l ectr i ca I activity of the heart from a first set of prefer- 
ably non-contact electrode sites and a second set of 

In general the present invention provides a method for in-contact electrode sites, 

producing a high-resolution, three-dimensional map of elec- jhe mapping catheter assembly 14 is coupled to interface 

trical activity of the inside surface of a heart chamber. 65 apparatus 22 which contains a signal generator 32, and 

The invention uses a specialized catheter system to obtain voltage acquisition apparatus 30. Preferably, in use, the 

the information necessary to generate such a map. signal generator 32 is used to measure the volumetric shape 
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of the heart chamber through impedance plethysmography. Catheter Description 

This signal generator is also used to determine the position FIG. 2 shows a portion of the mapping catheter assembly 

of the reference electrode within the heart chamber. Other 14 placed into a heart chamber 80. The mapping catheter 

techniques for characterizing the shape of the heart chamber assembly 14 includes a reference catheter 16 and an array 

may be substituted. Next, the signals from all the electrode 5 catheter 20. In FIG. 2 the array catheter 20 has been 

sites on the electrode array 19 are presented to the voltage expanded through the use of a stylet 92 to place the electrode 

acquisition apparatus 30 to derive a three-dimensional, array 19 into a stable and reproducible geometric shape. The 

instantaneous high resolution map of the electrical activity reference cat heter 16 has been passed through the lumen 82 

of the entire heart chamber volume. This map is calibrated of the calhetef 20 tQ lace a distal ti electrode 

by the use of a surface electrode 24. rhe calibration is both J(J ass ^ ^ ^ ^ an endocardial surface . In 

electrical and dimensional Lastly this three-dimensional ^ catheter des a mechanical location 

map, along with the signal from an intramural electrode 26 ' , . c \_ . A , in , 

preferably at the tip of the reference catheter 16, is used to refer ^ nce f f the P °?T" ° f lhe and ^ 

compute a two-dimensional map of the intramural electrical h P eleclrode assembly 24 provides an electrical potential 

activity within the heart wall. The two-dimensional map is reference at or in the heart wall for the mapping .process, 

a slice of the heart wall and represents the subsurface 15 Although the structures of FIG. 1 are preferred there are 

electrical activity in the heart wall itself. several alternatives within the scope of the invention. The 

Both of these "maps" can be followed over time which is principle objective of the preferred form of the catheter 

desirable. The true three-dimensional map also avoids the system is to reliably place a known collection of electrode 

problem of spatial averaging and generates an instantaneous, sites away from the endocardial surface, and one or more 

high resolution map of the electrical activity of the entire 20 electrode sites into contact with the endocardium. The array 

volume of the heart chamber and the endocardial surface. catheter is an illustrative structure for placing at least some 

This three-dimensional map is an order of magnitude more of the electrode sites away from the endocardial surface. The 

accurate and precise than previously obtained interpolation array catheter itself can be designed to mechanically posi- 

maps. The two-dimensional map of the intramural slice is tion one or more electrode sites on the endocardial surface, 

unavailable using prior techniques. 25 The reference catheter is a preferred structure for carrying 

Hardware Description one or more electrode sites and may be used to place these 

FIG. 1 shows the mapping system 10 coupled to a electrode sites into direct contact with the endocardial 

patient's heart 12. The mapping catheter assembly 14 is surface. 

inserted into a heart chamber and the reference electrode 24 It should be understood that the reference catheter could 

touches the endocardial surface 18. 30 be replaced with a fixed extension of the array catheter and 

The preferred array catheter 20 carries at least twenty-four used to push a segment of the array onto the endocardial 

individual electrode sites which are coupled to the interface surface. In this alternate embodiment the geometric shape of 

apparatus 22. The preferred reference catheter 16 is a coaxial the spherical array maintains the other electrodes out of 

extension of the array catheter 20. This reference catheter 16 contact with the endocardial surface, 

includes a surface electrode site 24 and a subsurface elec- 35 FIG. 3 shows the preferred construction of the mapping 

trode site 26 both of which are coupled to the interface catheter assembly 14 in exaggerated scale to clarify details 

apparatus 22. It should be understood that the electrode site of construction. In general, the array catheter 20 includes a 

24 can be located directly on the array catheter. The array flexible lead body 72 coupled to a deformable lead body 74. 

catheter 20 may be expanded into a known geometric shape, The deformable lead body 74 is preferably a braid 75 of 

preferably spherical Resolution is enhanced by the use of 40 insulated wires, several of which are shown as wire 93, wire 

larger sized spherical shapes. A balloon 77 or the like should 94, wire 95 and wire 96. An individual wire such as 93 may 

be incorporated under the electrode array 19 to exclude be traced in the figure from the electrical connection 79 at 

blood from the interior of the electrode array 19. The the proximal end 81 of the flexible lead body 72 through the 

spherical shape and exclusion of blood are not required for flexible lead body 72 to the distal braid ring 83 located on 

operability but they materially reduce the complexity of the 45 the deformable lead body 74. At a predetermined location in 

calculations required to generate the map displays. the deformable lead body 74 the insulation has been selec- 

The reference electrode 24 and/or the reference catheter tively removed from this wire 93 to form a representative 

16 serves several purposes. First they stabilize and maintain electrode site 84. Each of the several wires in the braid 75 

the array 19 at a known distance from a reference point on may potentially be used to form an electrode site. Preferably 

the endocardial surface 18 for calibration of the shape and 50 all of the typically twenty-four to one -hundred-twenty -eight 

volume calculations. Secondly, the surface electrode 24 is wires in the braid 75 are used to form electrode sites. Wires 

used to calibrate the electrical activity measurements of the not used as electrode sites provide mechanical support for 

endocardial surface 18 provided by the electrode array 19. the electrode array 19. In general, the electrode sites will be 

The interface apparatus 22 includes a switching assembly located equidistant from a center defined at the center of the 

28 which is a multiplexor to sequentially couple the various 55 spherical array. Other geometrical shapes are usable includ- 

electrode sites to the voltage acquisition apparatus 30, and ing ellipsoidal and the like. 

the signal generator apparatus 32. These devices are under The proximal end 81 of the mapping catheter assembly 14 

the control of a computer 34. The voltage acquisition has suitable electrical connection 79 for the individual wires 

apparatus 30 is preferably a 12 bit A to D converter. A signal connected to the various electrode sites. Similarly the proxi- 

generator 32 is also supplied to generate low current pulses 60 mal connector 79 can have a suitable electrical connection 

for determining the volume and shape of the endocardial for the distal tip electrode assembly 24 of the reference 

chamber using impedance plethysmography, and for deter- catheter 16 or the reference catheter 16 can use a separate 

mining the location of the reference catheter. connector. The distance 90 between the electrode array 19 

The computer 34 is preferably of the "workstation" class and the distal tip assembly 24 electrode can preferentially be 

to provide sufficient processing power to operate in essen- 65 varied by sliding the reference catheter through the lumen 

tially real time. This computer operates under the control of 82, as shown by motion arrow 85. This distance 90 may be 

software set forth in the flow chart of FIG. 12. "read" at the proximal end 81 by noting the relative position 
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of the end of the lead body 72 and the proximal end of the 
reference catheter 16. 

FIG. 4 is a view of the mapping catheter with the 
deformable lead body 74 in the collapsed position. 

FIG. 5 shows that the wire stylet 92 is attached to the 
distal braid ring 83 and positioned in the lumen 82. Traction 
applied to the distal braid ring 83 by relative motion of the 
stylet 92 with respect to the lead body 72 causes the braid 75 
to change shape. In general, traction causes the braid 75 to 
move from a generally cylindrical form seen in FIG. 4 to a 
generally spherical form seen best in FIG. 2 and FIG. 5. 

The preferred technique is to provide a stylet 92 which 
can be used to pull the braid 75 which will deploy the 
electrode array 19. However, other techniques may be used 
as well including an optional balloon 77 shown as in FIG. 3, 
which could be inflated under the electrode array 19 thereby 
causing the spherical deployment of the array 19. Modifi- 
cation of the braid 75 can be used to control the final shape 
of the array 19. For example an asymmetrical braid pattern 
using differing diameter wires within the braid can prefer- 
entially alter the shape of the array. The most important 
property of the geometric shape is that it spaces the electrode 
sites relatively far apart and that the shape be predictable 
with a high degree of accuracy. 

FIG. 6 shows a first embodiment of the reference catheter 
16 where the distal electrode assembly 24 is blunt and may 
be used to make a surface measurement against the endocar- 
dial surface. In this version of the catheter assembly the wire 
97 (FIG. 2) communicates to the distal tip electrode and this 
wire may be terminated in the connector 79. 

FIG. 8 shows an alternate reference catheter 98 which is 
preferred if both surface and/or subsurface measurements of 
the potential proximate the endocardial surface are desired. 
This catheter 98 includes both a reference electrode 24 and 
an extendable intramural electrode body 100. 

FIG. 9 illustrates the preferred use of an intramural 
electrode stylet 101 to retract the sharp intramural electrode 
body 100 into the reference catheter lead body 102. Motion 
of the intramural electrode body 100 into the lead body 102 
is shown by arrow 103. 

FIG. 10 shows the location of the intramural electrode site 
26 on the electrode body 100. It is desirable to use a 
relatively small electrode site to permit localization of the 
intramural electrical activity. 

The array catheter 20 may be made by any of a variety of 45 
techniques. In one method of manufacture, the braid 75 of 
insulated wires 93,94,95,96 can be encapsulated into a 
plastic material to form the flexible lead body 72. This 
plastic material can be any of various biocompatible com- 
pounds with polyurethane being preferred. The encapsula- 50 
tion material for the flexible lead body 72 is selected in part 
for its ability to be selectively removed to expose the 
insulated braid 75 to form the deformable lead body 74. The 
use of a braid 75 rather than a spiral wrap, axial wrap, or 
other configuration inherently strengthens and supports the 55 
electrodes due to the interlocking nature of the braid. This 
interlocking braid 75 also insures that, as the electrode array 

19 deploys, it does so with predictable dimensional control. 
This braid 75 structure also supports the array catheter 20 
and provides for the structural integrity of the array catheter 60 

20 where the encapsulating material has been removed. 
To form the deformable lead body 74 at the distal end of 

the array catheter 20, the encapsulating material can be 
removed by known techniques. In a preferred embodiment 
this removal is accomplished by mechanical removal of the 65 
encapsulating material by grinding or the like. It is also 
possible to remove the material with a solvent. If the 
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encapsulating material is polyurethane, tetrahydrofuran or 
cyclohexanone can be used as a solvent. In some embodi- 
ments the encapsulating material is not removed from the 
extreme distal tip to provide enhanced mechanical integrity 
forming a distal braid ring 83. 

With the insulated braid 75 exposed, to form the deform- 
able lead body 74 the electrodes sites can be formed by 
removing the insulation over the conductor in selected areas. 
Known techniques would involve mechanical, thermal or 
chemical removal of the insulation followed by identifica- 
tion of the appropriate conducting wire at the proximal 
connector 79. This method makes it difficult to have the 
orientation of the proximal conductors in a predictable 
repeatable manner. Color coding of the insulation to enable 
selection of the conductor/electrode is possible but is also 
difficult when large numbers of electrodes are required. 
Therefore it is preferred to select and form the electrode 
array through the use of high voltage electricity. By applying 
high voltage electricity (typically 1-3 KV) to the proximal 
end of the conductor and detecting this energy through the 
insulation it is possible to facilitate the creation of the 
electrode on a known conductor at a desired location. After 
localization, the electrode site can be created by removing 
insulation using standard means or by applying a higher 
voltage (eg. 5 KV) to break through the insulation. 

Modifications can be made to this mapping catheter 
assembly without departing from the teachings of the 
present invention. Accordingly the scope of the invention is 
only to be limited only by the accompanying claims. 
Software Description 

The illustrative method may be partitioned into nine steps 
as shown in FIG. 12. The partitioning of the step -wise 
sequence is done as an aid to explaining the invention and 
other equivalent partitioning can be readily substituted with- 
out departing from the scope of the invention. 

At step 41 the process begins. The illustrative process 
assumes that the electrode array assumes a known spherical 
shape within the heart chamber, and that there are at least 
twenty-four electrodes on the electrode array 19. This pre- 
ferred method can be readily modified to accommodate 
unknown and non-reproducible, non-spherical shaped 
arrays. The location of each of these electrode sites on the 
array surface is known from the mechanical configuration of 
the displayed array. A method of determining the location of 
the electrode array 19 and the location of the heart chamber 
walls (cardiac geometry) must be available. This geometry 
measurement (options include ultrasound or impedance 
plethysmography) is performed in step 41. If the reference 
catheter 16 is extended to the chamber wall 18 then its length 
can be used to calibrate the geometry measurements since 
the calculated distance can be compared to the reference 
catheter length. The geometry calculations are forced to 
converge on the known spacing represented by the physical 
dimensions of the catheters. In an alternative embodiment 
reference electrode 24 is positioned on array catheter 20 and 
therefore its position would be known. 

In step 42 the signals from all the electrode sites in the 
electrode array 19 are sampled by the A to D converter in the 
voltage acquisition apparatus 30. These measurements are 
stored in a digital file for later use in following steps. At this 
point (step 43) the known locations of all the electrodes on 
the electrode array 19 and the measured potentials at each 
electrode are used to create the intermediate parameters of 
the three-dimensional electrical activity map. This step uses 
field theory calculations presented in greater detail below. 
The components which are created in this step (0 /m ) are 
stored in a digital file for later use in following steps. 
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At the next stage the question is asked whether the 
reference catheter 16 is in a calibrating position. In the 
calibrating position, the reference catheter 16 projects 
directly out of the array catheter 20 establishing a length 
from the electrode array 19 which is a known distance from 
the wall 18 of the heart chamber. This calibration position 
may be confirmed using fluoroscopy. If the catheter is not in 
position then the process moves to step 45, 46 or 47. 

If the reference catheter 16 is in the calibrating position 
then in step 44 the exact position of the reference catheter 16 
is determined using the distance and orientation data from 
step 41. The available information includes position in space 
of the reference catheter 16 on the chamber wall 18 and the 
intermediate electrical activity map parameters of the three- 
dimensional map. Using these two sets of information the 
expected electrical activity at the reference catheter surface 
electrode site 24 is determined. The actual potential at this 
site 24 is measured from the reference catheter by the A to 
D converter in the voltage acquisition apparatus 30. Finally, 
a scale factor is adjusted which modifies the map calcula- 
tions to achieve calibrated results. This adjustment factor is 
used in all subsequent calculations of electrical activity. 

At step 47 the system polls the user to display a three- 
dimensional map. If such a map is desired then a method of 
displaying the electrical activity is first determined. Second 
an area, or volume is defined for which the electrical activity 
is to be viewed. Third a level of resolution is defined for this 
view of the electrical activity. Finally the electrical activity 
at all of the points defined by the display option, volume and 
resolution are computed using the field theory calculations 
and the adjustment factor mentioned above. These calcu- 
lated values are then used to display the data on computer 
34. 

FIG. 7 is a representative display 71 of the output of 
process 47. In the preferred presentation the heart is dis- 
played as a wire grid 36. The iso-potential map for example 
is overlaid on the wire grid 36 and several iso-potential lines 
such as iso-potential or isochrone line 38 are shown on the 
drawing. Typically the color of the wire grid 36 and the 
iso-potential or isochrone lines will be different to aid 
interpretation. The potentials may preferably be presented 
by a continuously filled color-scale rather than iso-potential 
or isochrone lines. The tightly closed iso-potential or iso- 
chrone line 39 may arise from an ectopic focus present at 
this location in the heart. In the representative display 71 of 45 
process 47 the mapping catheter assembly will not be 
shown. 

In step 45 a subthreshold pulse is supplied to the surface 
electrode 24 of the reference catheter 16 by the signal 
generator 32. In step 54 the voltages are measured at all of 50 
the electrode sites on the electrode array 19 by the voltage 
acquisition apparatus 30 One problem in locating the posi- 
tion of the subthreshold pulse is that other electrical activity 
may render it difficult to detect. To counteract this problem 
step 55 starts by subtracting the electrical activity which was 
just measured in step 44 from the measurements in step 54. 
The location of the tip of the reference catheter 16 (i.e. 
surface electrode 24), is found by first erforming the same 
field theory calculations of step 45 on this derived electrode 
data. Next, four positions in pace are defined which are 
positioned near the heart chamber walls. The potentials at 
these sites are calculated using the three-dimensional elec- 
trical activity map. These potentials are then used to 
triangulate, and thus determine, the position of the sub- 
threshold pulse at the surface electrode 24 of the reference 
catheter 16. If more accurate localization is desired then four 
more points which are much closer to the surface electrode 
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24 can be defined and the triangulation can be performed 
again. This procedure for locating the tip of the reference 
catheter 16 can be performed whether the surface electrode 
24 is touching the surface or is located in the blood volume 
and is not in contact with the endocardial surface. 

At step 48 the reference catheter's position in space can 
be displayed by superimposing it on the map of electrical 
activity created in step 47. An example of such a display 71 
is presented in FIG. 7. 

When step 46 is reached the surface electrode 24 is in a 
known position on the endocardial surface 18 of the heart 
chamber which is proper for determining the electrical 
activity of the tissue at that site. If the intramural or 
subsurface extension 100 which preferentially extends from 
the tip of the reference catheter 102 is not inserted into the 
tissue then the user of the system extends the subsurface 
electrode 26 into the wall 18. The potentials from the surface 
electrode 24 and from the intramural subsurface 26 electrode 
are measured by voltage acquisition apparatus 30. Next a 
line 21 along the heart chamber wall which has the surface 
electrode 24 at its center is defined by the user of the system. 
The three-dimensional map parameters from step 43 are then 
used to compute a umber of points along this line including 
the site of he reference catheter surface electrode 24. These 
alculations are adjusted to conform to the measured alue at 
the reference catheter surface electrode 24. ext a slice of 
tissue is defined and bounded by this line 21 (FIG. 7) and the 
location of the intramural subsurface electrode 26 (FIG. 11) 
and computed positions such as 23 and 25. Subsequently a 
two-dimensional map 27 of the electrical activity of this 
slice of tissue is computed using the center of gravity 
calculations detailed below in the section on algorithm 
descriptions. Points outside of the boundary of the slice 
cannot be computed accurately. In step 49 this map 27 of 
electrical activity within the two-dimensional slice is dis- 
played as illustrated in FIG, 11. In this instance the iso- 
potential line 17 indicates the location within the wall 18 of 
the ectopic focus. 

Description of the Preferred Computing Algorithms 

Two different algorithms are suitable for implementing 
different stages of the present invention. 

The algorithm used to derive the map of the electrical 
activity of the heart chamber employs electrostatic volume - 
conductor field theory to derive a high resolution map of the 
chamber volume. The second algorithm is able to estimate 
intramural electrical activity by interpolating between points 
on the endocardial surface and an intramural measurement 
using center of gravity calculations. 

In use, the preliminary process steps identify the position 
of the electrode array 19 consequently the field theory 
algorithm can be initialized with both contact and non- 
contact type data. This is one difference from the traditional 
prior art techniques which require either contact or non- 
contact for accurate results, but cannot accommodate both. 
This also permits the system to discern the difference 
between small regions of electrical activity close to the 
electrode array 19 from large regions of electrical activity 
further away from the electrode array 19. 

In the first algorithm, from electrostatic volume- 
conductor field theory it follows that all the electrodes 
within the solid angle view of every locus of electrical 
activity on the endocardial surface are integrated together to 
reconstruct the electrical activity at any given locus through- 
out the entire volume and upon the endocardium. Thus as 
best shown in FIG. 7 the signals from the electrode array 19 
on the catheter 20 produce a continuous map of the whole 
endocardium. This is another difference between the present 
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method and the traditional prior art approach which use the tial components anywhere in the cavity. Once again, the 
electrode with the lowest potential as the indicator of cardiac weighted components are summed to obtain the actual 
abnormality. By using the complete information in the potentials. Given that the potential is known on the probe 
algorithm, the resolution of the map shown in FIG. 7 is boundary, and given that the probe boundary is non- 
improved by at least a factor of ten over prior methods. 5 conductive, we can determine the coefficients A x and B lt 
Other improvements include: the ability to find the optimal yielding the following final solution for potential at any 
global minimum instead of sub-optimal local minima; the P oi ^ within the boundaries of the cavity, using a spherical 
elimination of blind spots between electrodes; the ability to P robe of radius R: 
detect abnormalities caused by multiple ectopic foci; the 

ability to distinguish between a localized focus of electrical 10 uV v[f— W ( 1 V-V"*" 1 "!* y (9 

activity at the endocardial surface and a distributed path of m , *°'Q£ I Ll2/ + lJlj?J + U/+lJW P" U r9) 

electrical activity in the more distant myocardium; and the 
ability to detect other types of electrical abnormalities 

including detection of ischemic or infarcted tissue. on exemplary method for evaluating the integral for 0 /BI 

The algorithm for creating the 3D map of the cardiac 15 is the technique of Filon integration with an estimating sum, 

volume takes advantage of the fact that myocardial electrical discretized by p latitudinal rows and q longitudinal columns 

activity instantaneously creates potential fields by electro- ^ electrodes on the spherical probe, 
tonic conduction. Since action potentials propagate several 

orders of magnitude slower than the speed of electrotonic An A A 

conduction, the potential field is quasi-static. Since there are 20 * ~jT q 4-{ 4^ Vi ° h ^ 

no significant charge sources in the blood volume, Laplace's ! 
Equation for potential completely describes the potential 

field in the blood volume: Note that p times q equals the total number of electrodes on 

the spherical probe array. The angle 8 ranges from zero to jc 

25 radians and ijj ranges from zero to 2x radians, 
^-o At this point the determination of the geometry of the 

endocardial walls enters into the algorithm. The potential of 

LaPlace's equation can be solved numerically or analyti- each poim on the en d OC ardial wall can now be computed by 

cally. Such numerical techniques include boundary element defining them as r, 0, and uj. During the activation sequence 

analysis and other interative approaches comprised of esti- 30 the gra phical representation of the electrical activity on the 

mating sums of nonlinear coefficients. endocardial surface can be slowed down by 30 to 40 times 

Specific analytical approaches can be developed based on t0 present a picture of the ventricular cavity within a time 

the shape of the probe (i.e. spherical, prolate spherical or frame useful for human viewing. 

cylindrical). From electrostatic field theory, the general A geometric description of the heart structure is required 

spherical harmonic series solution for potential is: 35 in order for the algorithm to account for the inherent effect 

of spatial averaging within the medium (blood). Spatial 

_yy ^ t _ u _ l} ' averaging is a function of both the conductive nature of the 

0(r, 9, <p) - 2^ 2j +Bir * v) medium as well as the physical dimensions of the medium. 

Given the above computed three-dimensional endocardial 

40 potential map, the intramural activation map of FIG. 11 is 

In spherical harmonics, Y /m (0, Uj) is the spherical bar- estimated by interpolating between the accurately computed 

monic series made up of Legendre Polynomials. <E /m is the endocardial potentials at locations 23 and 25 (FIG. 7), and 

lm** component of potential and is defined as: actual recorded endocardial value at the surface electrode 24 

and an actual recorded intramural value at the subsurface 

♦*,-Jf(0. WU*, n>)<*« 45 electrode 26 site. This first-order estimation of the myocar- 

where V(6, +) is the measured potential over the probe dial activation map assumes that the medium is homogenous 

radius R and dQ is the differential solid angle and, in and lh * * e medium contains no charge sources. This 

spherical coordinates, is defined as: myocardial activation estimation is limited by the fact that 

the myocardial medium is not homogeneous and that there 

dQ^s'm 6dd<h\> 50 are charge sources contained within the myocardial medium. 

If more than one intramural point was sampled the under- 

During the first step in the algorithmic determination of lying map of intramural electrical activity could be 

the 3D map of the electrical activity each <£lm component is improved by interpolating between the endocardial surface 

determined by integrating the potential at a given point with values and all the sample intramural values. The center-of- 

the spherical harmonic at that point with respect to the solid 55 gravity calculations can be summarized by the equation: 
angle element subtended from the origin to that point. This 
is an important aspect of the 3D map; its accuracy in creating « 
the 3D map is increased with increased numbers of elec- ^^017-73"*) 

trodes in the array and with increased size of the spherical v{7~ x ) = 

array. In practice it is necessary to compute the <£ /m com- 60 I l^-Tl"* 
ponents with the subscript set to 4 or greater. These 0 /m 4 
components are stored in an 1 by m array for later deter- 
mination of potentials anywhere in the volume within the where, VQ represents the potential at any desired point 
endocardial walls. defined by the three-dimensional vector x and, V, represents 
The bracketed expression of equation 1 (in terms of A lf 65 each of n known potentials at a point defined by the 
B lf and r) simply contains the extrapolation coefficients that three-dimensional vector t and, k is an exponent that matches 
weight the measured probe components to obtain the poten- the physical behavior of the tissue medium. 



